Probability measure generated by the superfidelity 
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We study the probability measure on the space of density matrices induced by the metric defined 
by using superfidelity. We give the formula for the probability density of eigenvalues. We also study 
some statistical properties of the set of density matrices equipped with the introduced measure and 
provide a method for generating density matrices according to the introduced measure. 

PACS numbers: 03.65.-w, 02.10.Yn, 45.10.Na 
Keywords: random states, quantum fidelity, superfidelity 



o 

1—5 



I 



> 

On 

(N 

l> 

O 



> 

X 



I. INTRODUCTION 

Recent applications of quantum mechanics are based 
on processing and transferring information encoded in 
quantum states. Random quantum states can be used 
to study various effects unique to quantum information 
theory pL, • This is especially true if one needs to get some 
information about the typical properties of the system in 
question [5]. In many cases it is important to quantify to 
what degree states are similar to the average state or how, 
on average, given quantity evolves during the execution 
of a quantum procedure. The crucial question emerging 
in this situation is how one should choose random sample 
from the set of quantum states. 

The aforementioned question be answered easily in the 
case of pure quantum states. In this situation there ex- 
ists a single, natural measure for constructing ensembles 
of states, namely Fubini-Study measure. The situation 
is more complex in the case of mixed quantum states. 
The probability measure can be introduced using vari- 
ous distance measures between quantum states [2 . By 
choosing the metric we also choose the probability mea- 
sure on the space of density matrices. Among the most 
commonly used metrics we can point out the trace dis- 
tance, Hilbert-Schmidt distance, and Bures distance. 

In the analysis of mixed quantum states Bures dis- 
tance is the most commonly used metric among the ones 
mentioned above. It has many important properties [2]. 
In particular it is a Riemannian and monotone metric. 
On the space of pure states it reduces to Fubini-Study 
metric [3] and it induces the statistical distance in the 
subspace of diagonal density matrices [? . 

The main aim of this paper is the analysis of the prob- 
ability measure on the space of density matrices induced 
by the metric defined in terms of superfidelity [5] . We cal- 
culate the formula for the probability density of eigenval- 
ues and study some properties of the space of quantum 
states equipped with the introduced measure. We also 
provide a method for sampling random density matrices 
according to the introduced distribution. 



This paper is organized as follows. In Section |TT] we 
introduce notation and basic facts used in the following 
sections. In Section ITTll we calculate the volume element 
for the measure generated by the metric based on superfi- 
delity and compare it with the analogous metric based on 



quantum fidelity. In Section IV we provide a formula for 
a probability density function on a simplex of eigenval- 
ues. We also calculate the normalization constant in the 
low-dimensional case. In Section Iv] we provide a method 
for sampling density matrices according to the introduced 
measure. Finally, in Section |VI| we provide a summary 
of the presented results. 



II. PRELIMINARIES 

Let use denote by A^jv the space of density matrices 
of size TV, i.e. TV X TV positive matrices with unit trace. 
By A we denote the simplex of eigenvalues. 

For two density matrices p,a G Mn, Bures distance 
can be defined in terms of quantum fidelity [3] as 



dB{p,'j)^j2-2^F{p,a), 



(1) 
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where fidelity, F{p,a) = [ti\y/py^\\ , provides a mea- 
sure of similarity on the space of density matrices. 

Probability measure on the simplex of eigenvalues gen- 
erated by Bures metric was calculated in [BH^. Various 
statistical properties of ensembles of quantum states with 
respect to this measure were discussed in \W . 

Bures distance is commonly used in quantum informa- 
tion theory as a natural metric on the space of density 
matrices. Unfortunately, fidelity used to express de has 
some serious drawbacks. In particular in order to cal- 
culate fidelity between two quantum states one needs to 
compute square root of matrix, which is in general com- 
putationally hard task. Also, fidelity cannot be measured 
directly in laboratory and thus cannot be used to analyse 
experiments directly. 

These drawbacks motivated the introduction of a new 
measure of similarity, namely superfidelity [5 , defined for 
p,cr & Mn as 



G(p, a) = tr per + v^l -tr/92VT 



trcr^ 



(2) 



Superfidelity shares many features with fidehty, i.e. it 
is bounded, symmetric and unitarly invariant. More- 
over it is jointly concave and supermultipHcative. It 
was proved that superfidehty gives an upper bound for 
fidehty, F{p,a) < G{p,a), where the equahty is for 
p,cr E M2 or in the case where one of the states is pure. 
It was also shown that, although G is not monotone [ID] . 
it can be used to define metric on A^jy. Using the cor- 
respondence between quantum operations and quantum 
states, superfidelity can be used to introduce metric on 
the space of quantum channels [TTj. Superfidelity was 
also proved to be useful in providing bounds on the trace 
distance jl2j {i.e. distinguisabiliy of states ^) and as a 
tool for studying new metrics on the space of quantum 
states [Ti] . 

In the following we use a metric on the space of density 
matrices defined for p,a G A4n as 



dG{p,a) = V^-2G{p,a). 



(3) 



Before we discuss further properties of this metric we 
should stress that the direct analogous of the Bures dis- 



tance, d /q(p, cr) = 1/2 — 2y/G{p, a), is not a metric. 

One should also note that, since G is not monotone it 
cannot be analysed using Morozova-Cencov-Petz theo- 
rem [2]. 



III. VOLUME ELEMENT FOR THE MEASURE 

To obtain the probability measure induced by metric 
Eq. (|3| one needs to derive the volume element. 

The calculations below follow the approach used by 
Hiibner [B]. We begin with the calculation of the line 
element 



4(p, p + dp) = 2- 2G{p, p + dp). 



(4) 



We introduce function A{t) = G{p, p + tdp), which allows 
to write the line element 



as 



1 d^ 
gijdp'dp' ^ 2 di2 [^g(p> P + * dp)] 



g,,dp'dp> = -A"{t) 



t=o 



(5) 



(6) 



Equivalently, with the use of matrix entries, the line ele- 
ment reads 

g,jdp dp> = i_y^2 + 2^\A{dp) \i). (7) 

Infinitesimal shift p + dp can be decomposed as a shift 
in eigenvalues and infinitesimal unitary rotation [7j 



p + dp^ p + dA+ [dU, p] , 



(8) 



where dA = '^dXi\i){i\ and {dUY = —dU. Rewriting 
dU in computational basis gives 



dU = ^{dxjk + idyjk)\j){k\ 



(9) 



with real coefficients dxjk — —dxkj and dyjk — dykj- 
After some calculations one gets 



trdp2 ^ J2idKf + 25](A, - Xjfiidx^.f + {dy,,f 



i<j 



and 



tr pdp = 2_, ^idXi. 



(10) 
(11) 



Expanding this we get the entries of the metric tensor 



A,:A,- 



g^jdp'dp^ = J2 ( ^ _\^ 2 + ^^i) d^^'dXj 

+2j2i^,^x,n{dx,,)^ + {dy.,n 



(12) 



1<J 



To obtain volume element of the sought measure, one 
must calculate the appropriate determinant 



dVa = Jdet ( j^^ + % ]dXi... dX„ 
X ]^2(Aj- Xjfdx.jdy, 

i<j 

Using the equality 



(13) 



XiXi 



ti-p^ 



1 



we obtain the expression for the volume element 
dAi . . . dA 



dVr,= 



yr^E;^ 



Y\_2{Xi- Xj)'^dxijdyij. (15) 



I I Kj 



One can compare the above formulas for the line el- 
ement with the analogous result for the metric given in 
terms of fidelity as 

dl,{p,p + dp) = 2{l-F{p,p + dp)). (16) 

In this case it is easy to check that the line element is 
given by formula 

dg,[p,p + dp) ^ V , , . ■ (17) 

'—^ Ai -f- Aj 

In the one-qubit case the above formula reads 
d\, {p, p + dp) 



2A(1-A) 



\dpiiV + \dp^2V + \dp2iV 

(18) 

where A and 1 — A are eigenvalues of p and dpij = {i\dp\j) 
and we have used the equality (l|dp|l) = — (2|dp|2). This 
is identical to Q for A^ = 2, which is what one expects 
since in this case F(p, a) = G{p, a). 



IV. PROBABILITY DENSITY FUNCTION 

In order to obtain probability measure we need to spec- 
ify the normalizing constant. This constant is an inverse 
of the integral of the volume element dVc over the group 
of unitary matrices and over the simplex of eigenvalues. 



over the simplex of eigenvalues. 

The above integral can be written in terms of expec- 
tation value with respect to Hilbert-Schmidt measure 



1 






E 



Vl-trp2 



(24) 



A. Normalization constant 

Integration with respect to U{N) is independent from 
the integration over the simplex of eigenvalues. We can 
rewrite Eq. [15] as 



where p is a random state distributed with Hilbert- 
Schmidt measure and 






r(iv2 



fN 



nLir(fc)r(fc + i) 



(25) 



dVr,^ 



The distribution of purity (trp^) for random states is 
r)N{N-i)/2 \ a subject of much study [15H12]- 

— - TT(^i — ^j)^ dXi . . .dXn\\ dxijdyij. The probability distribution function of purity is 



v/l-E.A,^ 



I i<j 



i^il 



After integrating this formula over U{N) we get 

r ( 2N{N-l)/2 



Vg = T 



N 



(19) 



^(A^-A,)Mc^Al...dA„, 



where Tjv is the volume of projective U{N) [H Eq. (148)] 

^N{N-l)/2 



-N 



Ul\'di 



(20) 
148)] 

(21) 



and A is the simplex of eigenvalues. 

Probability density function on a simplex of eigenval- 
ues is given by 



iGMW^C^YliK-X,)' 



1 



l<] 



vi-E,Ar 



(22) 



where Cn is a normalization constant. For A^ = 3 func- 



tion /cat is presented in Fig. 1(a) 



A A 



(a)Measure generated by 
\/l — G metric. 



(b)Measure generated by the 
Bures metric 



FIG. 1. Distribution of the eigenvalues for one-qutrit (A'^ — 3) 
density matrices for different probability measures. 

The normalization constant C^ is the following inte- 
gral 



h-'LW'"'''7r^: 



rG 



yr^E^A; 



:dX 



(23) 



known for Hilbert-Schmidt random states in the case of 
N — 2 and A^ = 3 IB]. Using these results, we can write 
explicitly normalizing constants 



C' 



2V2 



-C. 



HS 
2 ' 



G _ 432%/2 jjg 

^3 ~ ^Tt^T^^ • 



(26) 

(27) 



of 



In the case of A' > 3 one can use the series expansion 
2^ — and rewrite the above as 

\/l — r 



1 -^ oo 

'FiG ~ T^HS Z^ 



(2fc-l)!! 



2\fci 



E[tr(p^) 



(28) 



The moments of purity for Hilbert-Schmidt random state 
are given by [HI [T6] 



2\k^ 



E[(trp^) 



A^!(A2-1)! 



E 

felH \-k]s 



k\ 



V ' i.__l LI,^ =fe lii^l '*-!■ 



n 



{n + 2ki - i)l 
{q-i)W. 



Yl {2h-i-2kj+j). 



l<i<j<n 



(29) 



The constant C^ can be bounded from the above by 
using Jensen inequality 



1 
rG 



1 



E 



> 



1 

v/l-trp2 

1 



CW^l-E[trp^] CHS r 



2N 
N^ + 1 



thus 



HS 



'^N — L'Af 



1- 



2A^ 



(30) 
(31) 

(32) 



A2 + 1' 
Since the distribution of purity has the variance given by 

2fA2- 1)^ 
a2(trp2) = i^ 1 , (33) 

(A^2 + l)2^^2+2)(A2+3) 



it tends to be more concentrated around the mean given 

by 

2N 



E[trp^ 



N^ + V 



(34) 



which tends to zero for large N. For small x, function 
l/-\/l — X can be approximated with a small error by a 
linear function. Thus Jensen inequality gives good ap- 
proximation of C^ for large values of N, where tr p^ tends 
to be small. 



B. Mean purity 

Let pG be a random state distributed with measure G. 
Then the mean purity is given as 



Cat 






tip^ 



Vl-trp2 



(35) 



where p has Hilbert-Schmidt distribution. Next we have 



E 



tip^ 



Vl-trp2 



>E[trp2] E 



1 



Vl-trp2 



(36) 



which follows from the fact that random variables trp^ 
and , ^ are associated (see e.g. [18]). Finally, by 



using Eq. (24), we get 



C" 



E[tTpl]>-^E[trp']E 



v/l-trp2 



= E[trp2] 



, (37) 

From the above one can see that the mean purity for 
random state distributed with measure induced by the 
superfidelity is greater than the mean purity for random 
state distributed with Hilbert-Schmidt distribution. 



B. General case 

To generate random state of dimension N > 2 dis- 
tributed with measure induced by the superfidelity, one 
can use the rejection method (see e.g. [H]). 

Probability density function fcN on a simplex of 
eigenvalues can be bounded as 



/G.iv(A) < c/b,jv(A), VAgA 



(40) 



where fB,N is a probability density function generated by 
Bures measure [2] (see Fig. |l(b) ) 



/B,iv(A) — Cjy 

Indeed, we have 



VAi 



n 



{k-x,r 



. Xn \^ \i + Ao 



sup ■ 



/G,Ar(A)_CGiV-^/2(2/iV)^(^-l)/2 



and using the bound for C^ one can take 



(41) 



(42) 



0^T{N^)^^/^ 



nti r(i)2^(^-i)/2r(7V2/2)iv^'/2 



(43) 



as the constant in Eq. ( 40 1 



In order to generate a matrix distributed according 
to the measure induced by the superfidelity, one needs 
to generate a random matrix X distributed with Bures 
measure [20j and a random number u distributed uni- 
formly over the unit interval [0, 1]. To accept X as a ma- 
trix distributed according to the measure induced by the 
superfidelity, we check if w < ^ f "fx) holds. Unfortu- 
nately, constant c increases very rapidly with N and thus 
this method does not work very efficiently for large N . 



V. GENERATING RANDOM STATES 

A. One qubit case 

In the case of 2 x 2 matrices the density function on 
the simplex of eigenvalues reads 

2^2 



/g,2(A,1-A) = 



1 



TT v/A(l-A)' 



(38) 



Then the cumulative probability function for eigenval- 
ues by integrating /g,2 over interval [0, t] reads 

FgA^) = - (v/(l - i)t - 2^(1 - *)i^ + arcsin Vt) 

(39) 

From the above we obtain a simple method for gen- 
erating matrices with the above distribution. First one 
must generate eigenvalues of the matrix by inverting the 
cumulative distribution function and then rotate it by a 
random unitary matrix distributed with respect to Haar 
measure. 



VI. SUMMARY 

We have analysed random density matrices distributed 
according to probability measure induced by super- 
fildelity. We have derived the formula for the proba- 
bility density of eigenvalues according to this measure. 
We have also shown that random states distributed ac- 
cording to this measure have mean purity larger than in 
the case of Hilbert-Schmidt measure. We also provide a 
method for generating random matrices according to the 
introduced distribution. 

Still there are some problems which require further 
investigations. The first is the calculation of the exact 
formula for the normalization constant for the proba- 
bility density function. This is directly related to the 
distribution of purity for measures induced by the par- 
tial trace [15l [16] . The second problem is the inefficient 
method of sampling random states with the introduced 
measure, which could be used for numerical studies of 
the geometry of quantum states [5J [3TJ [^ . 
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